Multi-particle entanglement of hot trapped ions 
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We propose an efficient method to produce multi-particle entangled states of ions in an ion trap 
for which a wide range of interesting effects and applications have been suggested. Our preparation 
scheme exploits the collective vibrational motion of the ions, but it works in such a way that this 
motion need not be fully controlled in the experiment. The ions may, e.g., be in thermal motion 
and exchange mechanical energy with a surrounding heat bath without detrimental effects on the 
internal state preparation. Our scheme does not require access to the individual ions in the trap. 
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We describe in this Letter a method to construct states 
of the form 
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\y) = —{e^\gg...g)+e^\ee...e 



(1) 



where \gg...e) and |ee...e) (TV terms) are product states 
describing N ions which are all in the (same) internal 
state g or e. The unique correlation between different 
particles in such an entangled state is a significant non- 
classical feature of quantum mechanics. Entangled states 
with two particles have been employed to test Bell's in- 
equality and to rule out local realistic descriptions of Na- 
ture jlj . States like ([[]) with three particles, the so-called 
GHZ states Q, and with more particles allow a closer 
scrutiny of the role of quantum correlations, and there 
is a large current interest in carrying out experiments 
on such states. It has also been shown that the states 
(|I|) are ideal for spectroscopic investigations || , and the 
study of entangled states of many particles is in line with 
the rapidly growing interest in quantum encoding and 
quantum information processing Q. Together with an 
analysis of their applicability within high precision fre- 
quency measurements, a method to produce such states 
in an ion trap was proposed in || . That method requires 
experimental control of two collective vibrational modes 
of the ions, which are involved in stepwise manipulations 
of the internal states of the ions. 

In order to describe our more straightforward proposal, 
we shall show that the desired state (Q) can be obtained 
by application of one simple interaction Hamiltonian. We 
shall then show that this Hamiltonian is realized in an ion 
trap by simply illuminating all the ions with laser beams 
with two different colours, and we shall present numeri- 
cal results that display the performance of our proposal 
under quite relaxed assumptions for the functioning of 
the ion trap. 

We begin with all ions in the ground state \g), and we 
assume that our external perturbation interacts in the 
same way with all ions, so that the joint internal state re- 
mains symmetric with respect to interchange of the ions. 
In such a situation, it is convenient to apply the spin rep- 
resentation of atomic states, where a collective state of 
all ions is represented by | J, M) which is an eigenstate of 



the J z operator, J Z \J,M) = hM\J,M), M = - J,..., J, 
where h is Planck's constant. J z "counts" the number 
of excited ions, so that a state, invariant under permu- 
tation of the ions, with a certain number N e of excited 
ions is described as |J = N/2,M = N e - N/2), and 
the state (Q) is a superposition of the states \N/2, —N/2) 
and \N/2, N/2) . An increase of the quantum number M 
by unity corresponds to a single ionic excitation, and it 
is described by the angular momentum ladder operator 
J+\J,M) = hy/(J + M + 1)(J - M)\J, M + 1), where 
the normalization factor accounts precisely for the com- 
binatorial aspects of this excitation of the system. As we 
shall detail below, it is possible to apply classical laser ir- 
radiation which cannot be absorbed by a single ion, but 
which requires two ions to undergo the \g) «-> |e) transi- 
tion, simultaneously. When this field is applied to all ions 
in the trap, the effective Hamiltonian describing the evo- 
lution of the state of the ions in the spin representation 
can be written, 



J 2 _ + J+J- + J_J + ) 



4|4 



(2) 



where J_ is the hermitian conjugate of J + and where 
J x = (J + + J-)/2 denotes the x-component of the effec- 
tive spin. 

Time evolution is accommodated by application of the 
operator U(t) — exp(— i4% J 2 t/h 2 ) to the initial state 
| N/2, —N/2) and the population becomes distributed on 
all | J, M) states, with M differing from —N/2 by an even 
number. If N is even, population is transferred all the 
way to the state \N/2, N/2), which is the second com- 
ponent of (|l|). Moreover, at the instant t = ir/(8x), the 
state is precisely of the form ([!]), with <p g = — 7r/4 and 
(f> e = 7r/4 + Nir/2. We show in Fig. 1. the time evolution 
of the populations of the two extremal states for different 
(even) numbers of ions. At t = n/(8x) and at later odd 
multiples of this value, both populations equal one half. 

This result can be understood from the ro- 
tation properties of angular momenta: The ini- 
tial state can be expanded on eigenstates of J x : 
| N/2, -N/2) = J2m c m\N/2,M) x , and from the prop- 
erties of the Wigner rotation functions ||, it follows 
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that \N/2,N/2) = J2 M c M (-l) N/2 ~ M \N/2,M) x . In- 
sorting the values for <p g and <p e mentioned above, we 
can therefore write the state (|l|) as J2m c M-^(e _l7r / 4 + 
(— l) M e i7r / 4 )|iV/2, M) x . The net factors multiplying the 
initial amplitudes cm are unity for M even and — i for 
M odd. The action of U(t) in the J x basis amounts to a 
multiplication of each amplitude cj/ by exp(— 4ixM 2 t), 
and for t = tt/(8x) this factor just attains the desired 
value of unity for M even and —i for M odd. 

Note that the state (Q) is not only a multi-particle en- 
tangled state; it is a superposition of two "mesoscopi- 
cally distinguishable" states. In Quantum Optics, such 
states have been studied experimentally for single quan- 
tum systems, such as a harmonically trapped ion n and 
a single mode of the quantized radiation field j7j , excited 
into superpositions of states of "mesoscopic" separation. 
Mathematically, our analysis above resembles the one for 
production of such superpositions by propagation of an 
optical field through an amplitude dispersive medium . 

We now turn to a "microscopic" analysis of a collection 
of N two-level ions in a linear ion trap. The choice of a 
linear trap eliminates the so-called micro-motion for ions 
on the axis, and by choosing the transverse confinement 
much stronger than the longitudinal one, we can neglect 
coupling to transverse vibrations. Due to the Coulomb 
repulsion among the ions, the various longitudinal vibra- 
tional modes have different frequencies and by choosing 
the laser detuning close to the center-of-mass vibrational 
frequency we may exclude other modes from our analysis 

The quantized energy of the total system is composed 
of the center-of-mass vibrational energy of the ion string 
and the internal electronic energy of the ions. The vi- 
brational energy is represented by a harmonic oscillator 
with ladder operators a and and frequency v, and the 
internal energy is described by Pauli matrices a Z i and the 
energy difference Tiu) eg , 

H = hv{a) a + 1/2) + (3) 

i 

The ions interact with two laser fields with frequencies 
ojj (j = 1,2). We assume that both fields have the same 
intensity and we also assume that the ion-laser interac- 
tion strength characterized by the Rabi frequency Q is 
the same for all ions. Pauli matrices a±i describe inter- 
nal atomic transitions and the interaction Hamiltonian 
can be written 

H mt = ™(a +l e^^)S) + h .c). ( 4 ) 
i,3 

(h.c. in (^J) denotes the hermitian conjugate.) The Dicke 
parameters r)j are products of the projection of the k- 
vector of the j'th laser field along the direction of the 
string of ions and the rms excursion of the ionic center- 
of-mass along this direction jlOj . In practical realizations 
the A:-vectors can be nearly identical for the two fields, 



and we assume rji = r)2 = rj. Numerical solutions of 
the Schrodinger equation with the Hamiltonian given by 
Eqs.([|[i]), are readily obtained within the basis of 2 N 
internal states and a sufficiently large set of vibrational 
states {|/i/2---/jv)|")}, fi = e or g. First, however, we 
discuss a perturbative analysis pointing to appropriate 
choices of parameters which will effectively yield the spin 
Hamiltonian (g). 

The energy levels of the ion string constitute a lad- 
der of equidistant vibrational levels added to the internal 
ionic energies. This ladder structure has been extensively 
used in laser manipulation of ions in traps, where laser 
frequencies differing from the atomic transition frequency 
by a multiple of the trap frequency may be used to con- 
trol the internal and external dynamics H , and excitation 
of collective vibrational modes for several ions has been 
proposed as a means for implementing quantum logic in 
an ion trap quantum computer |Tl[| . 

Here, we make use of the idea we put forward in p5Q 
to apply laser fields with two different frequencies so that 
the two-photon process, exciting any pair of ions in the 
trap \gg) — » |ee) is resonant, i.e., uj\ + u>2 = 2uj eg , but 
neither of the frequencies are resonant with single ex- 
citations of the ions. By choosing io\ = uj eg — S and 
LU2 = uj eg + (5, where S is close to, but not resonant with, 
the center-of-mass vibration frequency we may limit our 
perturbative analysis to involve only the intermediate 
states with one excited ion and a vibrational quantum 
number raised or lowered by unity. Note that the lack of 
degeneracy of vibrational modes is crucial because it pre- 
vents a number of modes to interfere and obstruct the col- 
lective transitions in the multi-particle system. We show 
in Fig. 2 the effective transition paths between states 
\gg)\n) and |ee)|n) and between |eg)|n) and |ge)|n) me- 
diated by the interaction with the laser fields within any 
pair of ions. Although they serve as intermediate states 
in the process, the intermediate states are only virtually 
i.e., negligibly, populated. The interaction strengths of 
the different one-photon transitions in Fig. 2 depend on 
the value of the vibrational quantum number n and in the 
limit of r\\Jn + 1 << 1, terms between states n and n+1 
are proportional to \Jn + 1 whereas terms between n — 1 
and n are proportional to -Jn. When we add the con- 
tributions from the four different transition paths shown 
in either panel of Fig. 2, a "miracle" occurs: The paths 
(virtually) involving intermediate states with n+1 vibra- 
tion quanta yield [\Jn + l^lrj) 2 / (8 — v) and the ones with 
n — 1 quanta yield (^/nflrj) 2 / \v — 6). The signs of the de- 
nominators cause a destructive interference which exactly 
removes the n-dependence from the total effective Rabi 
frequency of the two-photon transition. This implies that 
the internal state evolves in a manner independent of the 
external state, which may be in a mixed (e.g., thermal) 
state and which may even evolve with time. 

When the summation over all ions is carried out, the 
couplings lead to the spin Hamiltonian (0) with 
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and c 2 = \/Tnthar ■ 
a mean excitation of nth 



and to the extent that our perturbative treatment is 
valid, we therefore expect to realize the results, presented 
in Fig. 1. In Fig. 3. we show numerical simulations of the 
dynamics governed by the correct Hamiltonian ([|]4|) for 4 
ions, i. e. beyond perturbation theory and to all orders in 
the Dicke parameter n. We assume a distribution on vi- 
brational levels with quantum numbers ranging between 
and 40, and we assume that the center-of-mass motion 
exchanges energy with a surrounding thermal reservoir 
represented by relaxation operators c\ = ^T(n t h + l)a 
Our calculation for a reservoir with 
5 vibrational quanta at the 
trap frequency and a coupling strength T = 0.0001 may 
represent a heating mechanism towards higher temper- 
atures counteracted by laser cooling on one ion in the 
string, reserved for this purpose. The action of the reser- 
voir is described by quantum jumps and a non-Hermitian 
Hamiltonian in a Monte Carlo wave function simulation 
fyl . Part (a) of the figure represents the result of a sin- 
gle Monte Carlo realisation with quantum jumps in the 
vibrational states at instants of time, indicated by arrows 
in the figure. Part (b) is an average over 10 such realisa- 
tions. The full lines present the populations and coher- 
ences of the internal state of the ions, obtained by tracing 
out the vibrational degrees of freedom. For comparison 
we also show the result of a calculation with the cor- 
responding spin Hamiltonian (^) (dashed curve). With 
the realistic parameters chosen, the state at vt ~ 1500 is 
very close to the desired state in (Q). Small oscillatory 
deviations from the smooth results of Fig. 1. are seen. 
The magnitude of the oscillations is given by the physical 
parameters, so that in the limit of small r\ and fi, they 
can be made arbitrarily small . 

The pairwise interaction with ions does not produce 
a coherent coupling of the two components in (|l|) if N 
is odd. In this case, however, one can verify within the 
spin model, that it suffices to apply a linear coupling 
Hi = 4:£J X for the duration t = 7r/(8£) in addition to 
application of our quadratic term ([2|). Since J x and J x 
commute, the linear Hamiltonian may actually be ap- 
plied before, after or simultaneously with the quadratic 
one. This additional operation is easy to implement in 
the ion trap, where one only needs to drive the ions on 
the atomic resonance frequency uj eg . 

In summary, we have presented a method to prepare a 
multi-particle entangled state of trapped ions. Our sug- 
gestion aims at experimental application, and it is rel- 
evant both for experiments with limited control of the 
vibrational motion, where the ions heat up for different 
technical reasons, or where there is not access to indi- 
vidual ions. The possibility to exercise control over some 
quantum degrees of freedom of a system while others are 
uncontrolled is common place in coherent optical excita- 
tion experiments, where the internal states of atoms in 
a gas may be controlled regardless of the external ther- 



mal motion. Our analysis adds the surprising possibil- 
ity of coherently controlling the mutual internal state of 
several particles by explicit use of the uncontrolled in- 
coherently excited degrees of freedom. We believe that 
much wider possibilities for quantum state control be- 
come feasible with the introduction of the crucial ele- 
ments of our method: resonance conditions allowing only 
pairwise transitions, virtual excitations of un-controlled 
degrees of freedom, and partially destructive interference 
between transition paths. We have recently applied these 
ideas to quantum computation in an ion trap |l2j and 
we are currently investigating their possible use in other 
physical systems. 
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FIG. 1. Time evolution of the population of the joint ionic 
ground state \gg---g) = l-^V/2, —N/2) (curves starting from 
the value of unity at t = 0), and the joint ionic excited state 
|ee...e) = \N/2, N/2) . Results are presented for different val- 
ues of the number of ions: N=2 (solid curves), N=4 (dashed 
curves) and N=8 (dotted curves). At t — n/(8x) the states 
are in a 50-50 superposition and the state (uh is obtained. 
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FIG. 2. Level scheme for a pair of ions sharing an oscillator 
degree of freedom. Left: by application of laser light with fre- 
quencies u>e S ± S, where S is somewhat smaller than the vibra- 
tional frequency v, we identify four transition paths between 
the states |<?<?)|n) and |ee)|n), which interfere as described in 
the text. Right: four similar transition paths are identified 
between states |e(?)|n) and |<?e)|n), yielding the same effective 
coupling among these states as between the states in the left 
panel. 




FIG. 3. Results of numerical integration of the Schrodinger 
equation for 4 ions. The physical parameters are 8 = 0.9v, 
D. = O.lv, rj = 0.1, T = O.OOOlu and n th = 5. (a) Single 
Monte Carlo realisation with quantum jumps in the vibra- 
tional states at instants of time, indicated by arrows in the 
figure, (b) Average over 10 realisations. The first full curve 
(counted from above at vt < 1000) is the population of the 
joint ground state, the second one is the imaginary part of 
the internal density matrix element between the joint ground 
and the joint excited state, the third one is the population of 
the joint excited state, and the fourth one is the real part of 
the off-diagonal density matrix element. The dashed curves 
are equivalent to the N = 4 curves of Fig. |l| obtained by 
application of the spin Hamiltonian (0) with \ — v ~ ~ 
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